Complex Square Well — A New Exactly Solvable Quantum Mechanical Model 
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Recently, a class of 'PT-invariant quantum mechanical models described by the non-Hermitian 
Hamiltonian H = p^ + x'^{ixy was studied. It was found that the energy levels for this theory are 
real for all e > 0. Here, the limit as e -+ oo is examined. It is shown that in this limit, the theory 
becomes exactly solvable. A generalization of this Hamiltonian, H = p^ +x'^^'^ (ix)" {M = 1, 2, 3, . . .) 
is also studied, and this PT-symmetric Hamiltonian becomes exactly solvable in the large-e limit as 
well. In effect, what is obtained in each case is a complex analog of the Hamiltonian for the square 
well potential. Expansions about the large-e limit are obtained. 
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I. INTRODUCTION 

The infinite square-well potential, 



the eigenvalues of H in Eq. (1.4) as functions of e are 
displayed in Fig. ^ 



^SW(2:) 



(|a;|<l), 

1 (1x1 = 1), 
oo (|a;| > 1), 



(1.1) 



is the simplest of all quantum potentials. It is studied 
at the beginning of any introductory class in quantum 
mechanics. This model is a useful teaching tool because 
the eigenvalues and eigenfunctions for this potential can 
all be found in closed form. 

The infinite square-well potential can be regarded as 
the limiting case of a class of potentials of the form 



Vm{x) = X 



2M 



(M= 1,2,3,4,...). 



Here, as M ^ oo, Vm{x) -^ Vsw(a;). 

The eigenvalues of the Hamiltonian for Vm , 



(1.2) 



(1.3) 



can only be found in closed form for the special case of 
the harmonic oscillator M = 1. For all other positive 
integer values of Af there is no exact solution to these 
anharmonic oscillators. Thus, the only two exactly solv- 
able cases known are the extreme lower and upper limits 
M — 1 and M -^ oo. The asymptotic behavior of the 
eigenvalues of H in Eq. (E3) for large M was studied in 
Ref. §. 

In a recent letter g] the spectra of the class of non- 
Hermitian 7^T-symmetric Hamiltonians of the form 



H = p^ + x^ (ix) 



(^>0) 



(1.4) 



were shown to be real and positive. It is believed that the 
reality and positivity of the spectra are a consequence of 
VT symmetry. Here, the case e = is again the har- 
monic oscillator. For finite values of e larger than there 
is no exact analytical solution for the eigenvalues. How- 
ever, solutions can be found by numerical integration; 




FIG. 1. Energy levels of the Hamiltonian _ff = p +x {ixY 
as functions of the parameter e. There are three regions: 
When e > the spectrum is entirely real and positive. All 
eigenvalues rise monotonically with increasing e. The lower 
bound of this region, e = 0, corresponds to the harmonic oscil- 
lator, whose energy levels are Ek = 2fc-|-l. When — 1 < e < 0, 
there are a finite number of real positive eigenvalues and an 
infinite number of complex conjugate pairs of eigenvalues. As 
e decreases from to —1, the number of real eigenvalues de- 
creases. As e approaches —1^, the ground-state energy di- 
verges. For e < — 1 there are no real eigenvalues. 

In the past we have always regarded the parameter e as 
being small; we have defined theories by analytically con- 
tinuing away from e = 0. However, in this paper w e in - 
vestigate the large-e limit of the Hamiltonian in Eq. (|l.4|). 
We will show that in this limit the theory becomes ex- 
actly solvable. An exact formula for the fcth energy level 
in the limit of large e is 



E.ie)^Uk+l 



(e^oo). 



(1.5) 



II. SPECIAL CASE M = 1. 



More generally, we will consider the large-e limit of an 
infinite number of classes of T^T-symmetric Hamiltonians 
of the form 

H=p^+ x^^'iixY (e > 0, M = 1, 2, 3, . . .)■ (1-6) 

For each positive integer value of A/, these Hamiltonians 
may be regarded as complex deformations of the Hermi- 
tian Hamiltonian H = p^ + x^*^ in Eq. (1.3). In the limit 
as e — > oo each of these Hamiltonians becomes exactly 
solvable; the spectrum for large e is given by 

Ek{M,e)^-^{k + jf^ e^ (e->cx)), (1.7) 

where P= 1,2,3,...,M. 

For the Hamiltonian H in Eq. (1.6) the Schrodinger dif- 
ferential equation corresponding to the eigenvalue prob- 
lem Hij} = E'4} is 



i}j"{x) + a;^*^ {ixyilj{x) = EiP{x). 



(1.8) 



To obtain real eigenvalues from this equation it is nec- 
essary to define the boundary conditions properly. The 
regions in the cut complex- a; plane in which ^{x) van- 
ishes exponentially as |a;| -^ cx) are wedges. In Refs. [|[^] 
the wedges for e > were chosen to be the analytic con- 
tinuations of the wedges for the anharmonic oscillator 
(e = 0), which are centered about the negative and posi- 
tive real axes and have angular opening -k /{M -f 1). This 
analytic continuation defines the boundary conditions in 
the complex- a; plane. For arbitrary e > the anti-Stokes' 
lines at the centers of the left and right wedges lie below 
the real axis at the angles 



eTT 



71cft = -TT - 



4A/-F2e-^4' 
eTT 



-'right 



4M -f 2e + 4 



(1.9) 



The opening angle of each of these wedges is 2-k/{2M + 
t + 2). In Refs. [||J3| the time- independent Schrodinger 
equation was integrated numerically inside the wedges 
to determine the eigenvalues to high precision. Observe 
that as e increases from its anharmonic oscillator value 
(e — 0), the wedges bounding the integration path un- 
dergo a continuous deformation as a function of e. As 
e increases, the opening angles of the wedges become 
smaller and both wedges rotate downward towards the 
negative-imaginary axis. Also, note that the angular dif- 
ference bright — ^loft = 27r(M -I- l)/e approaches zero as e 
increases. 

This paper is organized very simply. In Sec. ^ we 
consider the special case Af = 1 in Eq. (1.4). Then in 



Sec. [II we generalize to the case of arbitrary integer M 
in Eq. (1.6). Finally, in Sec. IV we examine expansions 
about the e ^ cx) limit of the theory. 



The eigenvalues of H in Eq. ( |l.4D can be found approx- 
imately using WKB theory. The left and right turning 
points for this calculation lie inside the left and right 
wedges at 



sicft = £;i/(2+^) exp -in 



,ght = i?^/(^+^)exp 



2e + 4 



2e + 4 



(2.1) 



As explained in Ref. p|, the WKB quantization formula 



k+ -\TT 



dx y/ E — x'^ {ixY (fc — > oo), 

(2.2) 



where the path of integration is a curve from the left 
turning point to the right turning point alo ng which the 
quantity dx times the integrand of Eq. ( ^.2[ ) is real. This 
path lies in the lower-half x plane and is symmetric with 
respect to the imaginary axis. The path resembles an 
inverted parabola; it emerges from the left turning point 
and rises monotonically until it crosses the imaginary 
axis; it then falls monotonically until it reaches the right 
turning point. As cal culated in Rcf. [E|, the WKB quan- 
tization formula (2^) gives 



Ek 






k 



i^Yi 



e+3 
e+2 



•■t+4 

€+4 



(fc-^oo). (2.3) 



W hen the parameter e is large, the right side of 
Eq. (2.3) simplifies dramatically and we have the result 
in Eq~(E5|) with corrections of order elne. As we will 
see, this happens to be the exact answer for all energy 
levels; that is, for all values of k. Since the WKB formula 



in Eq. (2^) is only valid for large k with e fixed, it is not 
at all obvious why the leading-order WKB calculation 
gives the exact answer. 

It is surprising to learn that the energy levels grow as 
e^ for large e. Recall that the energy levels of the Hamil- 



tonian in Eq. (1.3) approach finite limits as Af ^ oo. 
[These limits are the energy levels of the conventional 
square well Vsw in Eq. (^).] To understand why 
the energy levels for the PT-symmetric Hamiltonian in 
Eq. ( |1.4| ) grow as e^ we use the uncertainty principle. 
From Eq. ( |2.l| ) we see that the turning points rotate to- 
wards each other as e ^ cx). (They both approach the 
point —i on the negative imaginary axis.) Indeed, the 
distance between the turning points is o f or der 1/e. [For 
the case of the Hamiltonian H in Eq. ( pT^ ) the turning 
points stabilize at ±1 as M -^ cx).] Thus, the quantum 
particle is trapped in a region whose size Ax is of order 



1/e. The uncertainty in the momentum Ap of the parti- 
cle is therefore of order e. Finally, since the energy is the 
square of the momentum, we conclude that the energy 
levels must be of order e^. 

Let us rederive this result using the time-energy version 
of the miccrtainty principle. As explained in Refs. [§|,H, 
a classical particle described by the Hamiltonian in 



Eq. (1.4) exhibits periodic motion. The period T of this 
complex pendulum is given exactly by the formula 



T = \^E- 



^r(f±f)cos(3fL_ 

U+2e 



For large e we have 



T - 47r/(eV£;) (e -» c»). 



(2.4) 



(2.5) 



Multiplying this equation by E gives the product ET on 
the left side, which, by the uncertainty principle is of 
order 1. Thus, solving for E^ we find again that E is of 
order e^ for large e. 

This last calculation illustrates an important differ- 
ence between conventional quantum theories and VT- 
symmetric quantum theories. In a conventional Hermi- 
tian theory both the classical periodic motion and the 
WKB path of integration coincide; this classically al- 
lowed region lies on the real axis between the turning 
points. For 'PT-symmetric theories the WKB contour 
and the classical path do not coincide. The classical pe- 
riodic motion follows a path joining the turning points 
that, like the WKB path, is symmetric about the nega- 
tive imaginary axis. However, unlike the WKB path, the 
classical path moves downward rather than upward as it 
approaches the negative-imaginary axis (see, for exam- 
ple, Fig. 2ofRef. §). 

Having discussed this problem hcuristically, we now 
give a precise calculation of the spectrum in the limit of 
large e. We begin by substituting 



£'5+7 



(2.6) 



into Eq. (L8) with Af = 1. The resulting differential 
equation for large e is 



—^^{z) + F7^^1 + e^^^)i;{z)^0, 



where 



F = E/e^ 



(2.7) 



(2. 



and we have used the identity lime^oo(l + x/eY — e^. 

The advantage of the differential equation ( ^.7| ) is that 
it is independent of e. [As such, this equation corre- 
sponds to the M-independent Schrodingcr equation for 
the square well that is obtained from Eq. (1.3) in the 
limit of large A/.] In the variable z the turning points 



at z = — 1 and z = 1 are fixed and well separated in the 
limit of large e. The larg e-e b ehavior of E in Eq. (|l.5| ) 



is already evident in Eq. (2.8). Imposing the appropri 



ate boundary conditions on Eq. (2/7) gives eigenvalues F 
that are clearly independent of e. Thus, for large e we 
see that E grows like e^. 

Be cau se there is no longer any small parameter in 
Eq. (2.7), this equation cannot be solved approximately 
using a perturbative method such as WKB. It is neces- 
sary to solve this equation exactly. Fortunately, we can 
solve it exactly by making a simple substitution. The 
change of variable 



= 2\/Fe"^/2 



(2.9) 



converts Eq. ( |2.7| ) to a modified Bessel equation [^1 

w'^-—r^l;(w) + w-—^p{w) - iv? + i/^) ip{w) = 0, (2.10) 
dw^ dw 



where 



v = 2%/F. 



(2.11) 



The exact solution to this equation is a linear combina- 
tion of modified Bessel functions ptl : 



^[w) = Cil,{w) + C2K,{w), 



(2.12) 



where Ci and C2 are arbitrary constants. Thus, in terms 
of the z variable we have 

^{z) = Cil, (ve"'-'"^^ + C2K, (ye'^'/^^ . (2.13) 

We must now impose boundary conditions on '4){z). 
Emanating from the turning points at z = — 1 and z = 1 
are three Stokes' lines (lines along which the solution is 
purely oscillatory and not growing or falling exponen- 
tially) and three anti-Stokes' lines (lines along which the 
solution is purely exponential and not oscillatory) . These 
Stokes' and anti-Stokes' lines are shown as dashed and 
solid lines on Fig. g. The Stokes' lines emerge from the 
turning points going up to the left and the right at 30° 
and also directly down. The anti-Stokes' lines emerge 
from the turning points going down to the left and the 
right at 30° and directly up. Note that the Stokes' line 
going up to the right from the turning point at z = — 1 
joins continuously onto the Stokes' line going up to the 
left from the turning point at z = 1. The anti-Stokes' 
lines going down to the left from z = — 1 and down to the 
right from z = 1 eventually become vertical and asymp- 
tote to the lines Re z = — 2 and Re z = 2. We impose the 
boundary conditions that '4){z) ^ on these anti-Stokes' 
lines because these correspond to the center lines of the 
wedges in Eq. (|l.9| ) in the complex- a; plane (for M = 1). 

To s umm arize, in the large-e limit of the Hamiltonian 
in Eq. ( |l.4| ), the eigenvalue problem for the scaled eigen- 
values _F is a two-turning-point problem that lies along 
an arch-shaped contour. The legs of the arches lie below 



the real-z axis and approach ±2— ioo. The turning points 
at z = ±1 are joined by the Stokes' hne lying above the 
real-z axis as indicated in Fig. g. This is the complex 
version of the infinite-square-well problem in elementary 
quantum mechanics. In the square-well problem there 
are also two turning points at ±1 joined by a Stokes' line 
lying on the real axis. However, there are no anti-Stokes' 
lines along which the wave function dies away exponen- 
tially; the wave function simply vanishes at the turning 
points. 

The quantized energy levels are determined by impos- 
ing the boundary conditions discussed above on the mod- 
ified Bessel functions in Eq. (2.13). For simplicity, we 
impose these conditions on the vertical lines z = ±2 — ly, 
where y — > -|-cx3. In te rms of the variable y the wave 
function tp in Eq. (2.13) becomes 



4,{y) = Ci/, ( 



j,g-*^g^y/2 



C.K, 



j^g-i^g^2//2 



(2.14) 



at z = —2 — iy, and 

^{y) = Cih (ye'^'e'^y/^) + C2K, (ve'^e^y'"^^ (2.15) 
a,t z — 2 — iy. 



FIG 




Stokes' lines and anti-Stokes' lines for the differ- 



ential equation (2.7). Three Stokes' lines (dashed lines) and 
three anti-Stokes' lines (solid lines) emerge from the turning 
points at 2 = ±1. The pa th o f integration for the WKB quan- 
tization condition in Eq. (2.2) corresponds to the arch-shaped 



dotted hne connecting the turning points. 

Our objective now is to simplify these equations by 
making the arguments of the modified Bessel equations 
entirely real and positive. To do so we use the following 
functional equations satisfied by I^, and Kn [Q] : 

J^(e™"2) =e™'''^*/^(z), 

K, (e™"z) = e~"^''^^K,{z) ~ z^!^^i^/,(z), (2.16) 

sin(i/7r) 

wher e m is an integer. According to these relations, 
Eq. (2.14) becomes 



^{y) = Cie--^"/, (ve^y/^) 



+ C2 



e'"''K, 



^ (ve'^y'A + mh ive'^y'A] (2.17) 



and Eq. ( 2.15 ) becomes 



V;(y) = Cie'''^V. (i/e-^/2 



C, 



K^ (ve-^y/A - ml^ [ve-^y/A] . (2.18) 



Next, we use the asymptotic behavior of the modified 
Bessel functions for large positive argument. The func- 
tion /i/(r) grows exponentially and the function K^ir) 
decays exponentially for large positive r H| : 

Iy (r) ^ e** (r->+oo), 




(r — > -l-oo). 



(2.19) 



Eliminating the growing exponentials in Eqs. (2.17) and 
( ^.18 ) gives a pair of linear equations to be satisfied by 
the coefficients Ci and C2: 



Cie-'^'^^ + C2iTT = 0, 



(2.20) 

A nontrivial solution to Eq. ( 2.2C ) exists only if the 
determinant of the coefficients vanishes: 



det 



-2z7rcos(z^7r) — 0. 



(2.21) 



TABLE I. Compari son o f the numerical values of F{e) with 
that predicted in Eq. (2.22) for the ground state fc = of an 
x^{ixy theory. The second column gives the exact values of 
the ground-state energy for various values of e in the first col- 
umn. In the third column is the value of F obtained from the 
exact energy in the second column using Eq. (Hdl), which is a 
more precise version of Eq. {2.i). Finally, in the fourth and 
fifth columns are the first and second Richardson extrapolants 
P of the numbers in the third column. Note that the exact 
values of F{e) and their Richardson extrapolants rapidly ap- 
proach the asymptotic value 1/16 = 0.0625. 



Eoje) 



F{e) 



Ri{^) 



M^) 



8 5.55331 




0.07825 


- 






- 


18 20.67629 




0.06998 


0.06336 


- 


28 46.94324 




0.06742 


0.06281 


0.06259 


38 84.78728 




0.06617 


0.06266 


0.06253 


48 134.43752 




0.06542 


0.06260 


0.06251 


58 196.03417 




0.06493 


0.06257 


0.06251 


Hence, v — fc -1- ^, 


and from 


Eq. ( 


2.11 


), we 


have the 


exact result 










F=-{k + 
4 V 


^ 


(fc- 


0,1,2 


3,. 


..). 


(2.22) 



Finally, we use Eq. ( |2.8| ) to obtain the large-e behavior 
of the eigenvalues E given in Eq. (1.5). We verify this 
result numerically in Table |. 



III. ARBITRARY INTEGER M 

The calculation of the ener gy levels for the general class 
of theories given in Eq. (1.6) is a straightforward gener- 
alization of the calculation for the case M = 1 in Sec. ||. 
The crucial ingredient in the calculation is understanding 
the array of Stokes' and anti-Stokes' lines along which we 
impose the boundary conditions. This difference leads to 
M-dependent wave functions, but the condition that de- 
termines the eigenvalues is still a simple trigonometric 
equation. 

Just as for the case M = 1, we scale the differential 
equation (l.S) using Eq. (|2.6|). In the limit as e — > cx 3 



the resulting differential equation is identical to Eq. (2.7) 
except that now there is a factor of (—1)^^+^ multiplying 
the exponential term. Again, we define F as in Eq. (|2.S| ) 
and change to the variable w as prescribed by Eq. (2.9). 
This gives the differential equation 



w 



dw^ 






w——%Ij{w) 
aw 



'w^ 



+ v^] iIj{w) = 0, 



(3.1) 

which is the generalization of Eq. (^.IQ ). In this equation 
V is defined as before by Eq. ( [2.11 ) . Note that except for 



the appearance of the (— 1)*^+^ multiplying the vo^ term 
this equation is independent of M . 

To solve Eq. (|]j) we consider the two cases of odd 
M and even M s epar ately. If M is odd this equation is 
identical to Eq. (|2.10|), and the general solution is that 



given in Eq. (2.12). If M is even Eq. (3.1) is no longer a 



modified Bessel equation, but instead is just the standard 
Bessel equa tion. Hence, in this case the general solution 
to Eq. ( [3.l[ ) is a linear combination of the ordinary Bessel 
functions Ji, and Y^: 



i:{w) = CiJ,{w) + C2Y,{w). 



(3.2) 



Thus, in terms of the variable z the wave function in this 
case is 

i^{z) = CiJ, (j/e""/2) + C2Y, (i/e*""/2^ . (3.3) 

Although it appears that this solution is independent 
of the parameter M , one must recall that the boundary 
conditions do depend on M. Thus, the wave functions 
and energy eigenvalues do indeed depend on M . To be 
precise, for a given M the Stokes' lines emanating from 
z = ±M are joined by a string of adjacent arches of 
length 2. The anti-Stokes' lines leave ±M and asymptote 
to the lines Rez = ±(M + 1). 

For odd M we impose the boundary conditions as for 
the case M = 1 except that the wave function -0 vanishes 
along different lines. For even M the wave functions in 
Eq. ( |3.3| ) must first be expressed in terms of modified 
Bessel functions using the functional equations B 



Jviiz) 



'/'Liz 



Y,{iz) = —e-''^''^K,{z) + ie-^''''h{z), 



(3.4) 



and then be treated using the same procedure as for odd 
M. 

Although the matrix elements for the linear equations 
obtained for odd M and even M are quite different, the 
eigenvalue conditions are similar. For M = 2 the condi- 
tion is 



cos(2j/7r) = -1/2, 



whose solution is 



and 



. = k+-. 



Thus, for large e the energy is given by 

1 / P^' 
^^4^+3 



(3.5) 



(3.6) 



(3.7) 



whe re P =1,2. Note that this result is the M = 2 case of 
Eq. ( |1.7D . This expression is verified numerically in Table 
y for the case of the ground state energy corresponding 
to A; = and P — 1. For arbitrary M one obtains the 
result in Eq. (1.7). 



TABLE II. Comparison of t he n umerical values of P(e) 
with that predicted in Eq. (3.7) for the ground state 



fc = 0, P = 1 of an a; {ixY theory. The second column gives 
the exact values of the ground-state energy for various val- 
ues of e in the first column. In the third column is the value 
of F obtained from the energy in the second column using 
Eq. (kj). Finally, in the fourth and fifth columns are the 
first and second Richardson extrapolants [H] of the numbers 
in the third column. Note that the exact values of P(e) and 
their Richardson extrapolants rapidly approach the asymp- 
totic value 1/36 = 0.0277778. 



e 


^o(e) 


P(e) 


Pi (6) 


R2[e) 


8 


2.65128 


0.05035 


- 


- 


18 


9.21477 


0.03551 


0.02661 


- 


28 


20.70525 


0.03232 


0.02722 


0.02740 


38 


37.32010 


0.03097 


0.02746 


0.02766 


48 


59.16865 


0.03023 


0.02756 


0.02772 


58 


86.31766 


0.02977 


0.02764 


0.02775 



Observe that the magnitude of the energy eigenvalues 
decreases as M increases. At first glance this might seem 
surprising, but it can be easily understood in terms of the 
uncertainty principle. As M increases, the anti-Stokes' 
lines on which we impos e th e boundary conditions for 
the differential equation (1.8) move away fro m th e neg- 
ative imaginary axis, as we can see from Eq. (l.£). For 
example, for fixed e the anti-Stokes' lines for M ~ 2 are 
separated by a greater distance than for M = 1; the 
anti-Stoke's lines for M — 3 are separated by a greater 
distance than for M = 2, and so on. Hence, the un- 
certainty in the position Ax increases with M. By the 
uncertainty principle, this increase in the uncertainty of 
the position corresponds to a decrease in the uncertainty 
of the momentum, and thus, a decrease in the energy. 
This argument explains the large-M behavior of the re- 
suh in Eq. 1^. 



IV. HIGHER-ORDER CORRECTIONS TO THE 
t -> oo LIMIT FOR M = 1 

In this section we show how to calculate the corrections 
to the large-e behavior in Eq. (f^). These corrections are 
of order e and e In e. From these higher-order calculations 
we obtain an extremely accurate approximation Ek for 
all k. Our asym ptot ic analysis begins with the change of 
variabl e in Eq. ( 2^ ) , but we use a more precise version 



F{e) 



€ + 4 

E^+^ 



(4.1) 



We find that the function F{e) is a series in inverse pow- 
ers of e of the form F = fa + fie~^ + /ae'^ + ■■■. The 



coefficient /o — F{oo) is given in Eq. (2.22). Our objec- 
tive here is to calculate /i , and from this to calculate the 
first correction to E. 

In addition to F{e), the wave function ip is also a series 
in inverse powers of e, il>{z) = 'ipo + 'iJji{z)e~^ + il^2{z)e~'^ + 
■ ■ ■. Using this series and collecting like powers of e we 
obtain the following sequence of differential equations: 

. _^ 
■ dz2' 

fP 
e-i : — , 

dz2 



^ipo{z) + /o7r^ (l 



1 + e' 



Vo(2) = 0, 



/„^2g-- +/,(! 



i^M^) + foT:' {1 + e^'^) Mz) 



dz 






he' 



+fie'- 



,2„2 



./l(l 



3„3 ■ 



Tr'^Tpoiz), 



n'^ipiiz) 



z TT + z^n-^i 



2^2 



.zw 



+ ,f2 (1 + e'^'^) 



TT^Mz)- (4.2) 



The first equation is exactly Eq. (2.7). The second 
equation contains the coefficient /i. To solve for /i we 
observe that the solution to the homogeneous part of 
the equation is just the solution to the first equation. 
This suggests using the method of reduction of order; to 
wit, we let tjji{z) = Ui{z)'ipo{z) ■ To solve the resulting 
equation for /i we then multiply by ipoiz) and integrate 
over the WKB path with respect to z. The first equation 
can be used to simplify this equation and the left side 
becomes the expression ui(z)V'o(2) evaluated at the end 
points. Since ipQ^—ioo) = 0, the left side equals zero, and 
we can solve for /i in quadrature form. To be explicit. 



_ 1 , /!— dzzVg(z)e- 

/I - ^2-'°'^ r2-ioo 



/: 



2 — 100 



dzV^'(^)(l 



(4.3) 



To prepare for evalu atin g these integrals we change to 
the variable w in Eq. (^.9[) with F = fo and obtain 



/i = 



i/- 



-CX3 — 2(5 

oo-\-i5 



dw wipQ (w) In 



L 



oo-\~ iS 



1f^o^H(i + ^) 



(4.4) 



where S is infinitesimal and the contour of integration 
goes around the origin. 

For the case fc = these integrals are easy to evaluate 
because /o = 1/16 and i/'o(u') = Ii/2{w) + Ki/2{w)/it 



e^/v^TTu;. Substituting these expressions into Eq. ( |4.4| ) 
gives 



1 r°^-''dwe'^\n\2w) 

p ^ J— OO + ld ^ '' 

/1-2 r-oo-.5^ 2»(4+J^) 



(4.5) 



By carefully evaluating the discontinuities across the 
branch cut and the residues at the singularities of the 
integrands, we obtain 



/i = 7/4, 



(4.6) 



wher e 7 is Euler's constant. Combining this result with 
Eq. ( O ) and solving for E in the limit of large e yields 



E= — e2-ielne+i(l + 7 + 21n2)e + 0(lne). (4.7) 
16 4 4 

By comparison, if we calculate to next order in WKB, 
we obtain for the kt\i energy level 



Ek 



r(!^)V^(fc + i/2) 



2+e 



m 



4-f-2e 
4+e 



(2 + e)(l 



^in ( JT^) 



67r(fc + 1/2)2(4 -he)2 



(fc -^ 00). (4., 



Taking the large e limit of this expression gives 

^=^-^flne+iQ+ln2^e + 0(lne). (4.9) 

The appearance of a log e term in this behavior is a con- 
sequence of the structure of Eq. (4.1). Note that the 
coefficent of the e term for WKB differs from the ex- 
act result but WKB is numerically very accurate. WKB 
gives 0.75662 compared with 0.74088 for the exact result. 
In Table III the results of a Richardson extrapolation 



of the exact values of F{e) are given. These results 
verify that the value of /i is correct. 
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TABLE IIL Co mpa rison of the numerical value of the co- 



efficient /i in Eq. (4.6) with a fit to the exact values of F{e) 
for the case of the ground state fc = of an x^{ixy theory. 
The second column gives the exact values of the ground-state 
energy for various values of e in the first column. In the third 
column is the approximation to /i obtained from F(e)_by sub- 



tracting off the leading large-e behavior given in Eq. (2.22). In 
the fourth and fifth columns are the first and second Richard- 
son extrapolants fel of the numbers in the third column. 
Note that the approximations in column 3 and their Richard- 
son extrapolants rapidly approach the asymptotic value of 
/i =7/4 = 0.144304. 



e 


Eo{e) 


Roie) 


Ri{e) 


R2{e) 


8 


5.55331 


0.12597 


- 


- 


18 


20.67629 


0.13460 


0.14150 


- 


28 


46.94324 


0.13767 


0.14321 


0.14389 


38 


84.78728 


0.13926 


0.14372 


0.14418 


48 


134.43752 


0.14024 


0.14394 


0.14425 


58 


196.03417 


0.14090 


0.14406 


0.14428 



